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ABSTRACT 

We present the results of a dedicated campaign on the afterglow of 
GRB 030329 with the millimeter interferometers of the Owens Valley Radio Ob- 
servatory (OVRO), the Berkeley-IUinois-Maryland Association (BIMA), and with 
the MAMBO-2 bolometer array on the IRAM 30-m telescope. These observations 
allow us to trace the full evolution of the afterglow of GRB 030329 at frequencies 
of 100 GHz and 250 GHz for the first time. The millimeter light curves exhibit 
two main features: a bright, constant flux density portion and a steep power- 
law decline. The absence of bright, short-lived millimeter emission is used to 
show that the GRB central engine was not actively injecting energy well after 
the burst. The millimeter data support a model, advocated by Berger et al., of a 
two-component jet-like outflow in which a narrow angle jet is responsible for the 
high energy emission and early optical afterglow, and a wide-angle jet carrying 
most of the energy is powering the radio and late optical afterglow emission. 

Subject headings: gamma rays: bursts — radio continuum: general — cosmology: 
observations 



1. Introduction 

A link between long-duration gamma-ray bursts (GRBs) and the core collapse of mas- 
sive stars has long been claimed on observational (e.g. Bloom, Kulkarni & Djorgovski 2002; 
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Bloom et al. 2002; Price et al. 2002) and theoretical grounds (Woosley 1993; Paczynski 1998; 
MacFadyen & Woosley 1999), but the recent GRB 030329 has strengthened this association 
considerably. Optical spectra taken of this event (Stanek et al. 2003; Hjorth et al. 2003; 
Kawabata et al. 2003) showed the usual power-law continuum from the afterglow, super- 
imposed upon which were lines characteristic of a Type Ic supernova (SN). Designated 
SN 2003dh, the brightness of this SNe and its broad line widths are strikingly similar to 
another pecuhar Ic SN 1998bw (Patat et al. 2001; Maeda et al. 2002), or perhaps SN1997ef 
(Iwamoto et al. 2000). Depending on the degree of asphericity assumed in modehng the 
explosion, the derived kinetic energies for SN 1998bw and SN 1997ef are in the range of 
5-50 foe (1 foe=10^^ erg). Such events are labeled "hypernovae" , an empirical classifi- 
cation to distinguish them from ordinary core collapse SNe with energies of order 1 foe 
(Iwamoto et al. 1998). 

While the non-relativistic component of the explosion (as traced by SN 2003dh) may 
have been hyper-energetic, the relativistic component (as traced by GRB 030329) seems to 
be sub-energetic. As noted by Granot, Nakar & Piran (2003) and Berger et al. (2003), if the 
sharp break in the optical light curves at t = 0.55 days (Price et al. 2003) is due to a jet-like 
outflow, then the gamma-ray energy released ~0.05 foe and the X-ray luminosity (at t=10 
hrs) Lx ~ 3 X 10"^^ erg s~^ . These values of and Lx lie an order of magnitude or more below 
which most GRBs are tightly clustered (Bloom, Prail & Kulkarni 2003; Berger, Kulkarni & 
Frail 2003). Outliers in the energy/luminosity distribution are potentially important for 
exploring the diversity of the GRB phenomena, and how the central engine partitions the 
explosion energy. 

Millimeter detections of gamma ray bursts (GRBs), while difficult to achieve (e.g., 
Bremer et al. 1998, Shepherd et al. 1998), are a potentially powerful diagnostic of the 
explosion energy. Since the peak of the synchrotron spectrum is expected to pass through 
the millimeter band on a time scale of a day or so after the burst (Sari, Piran & Narayan 
1998), such observations are useful for constraining the peak of the spectrum, a quantity that 
is difficult to obtain by other means. When combined with broadband data these millimeter 
observations have been especially useful in deriving the kinetic energy of the outflow and the 
density structure of the circumburst environment (Galama et al. 2000; Berger et al. 2000; 
Harrison et al. 2001; Yost et al. 2002). 

In this paper we present measurements of GRB 030329 at 100 GHz made with the 
millimeter interferometers of the Owens Valley Radio Observatory (OVRO) and the Bcrkeley- 
lUinois-Maryland Association (BIMA), and measurements at 250 GHz made with MAMBO-2 
bolometer array on the IRAM 30 m telescope. At 2;=0.1686 (Caldwell et al. 2003) this is the 
closest known GRB and subsequently the flux density at millimeter wavelengths was more 
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than ten times brighter than any previous event. These observations allow us to trace the 
full evolution of the afterglow at millimeter wavelengths for the first time. In §2 we discuss 
the observations and calibration issues, in §3 and §4 we present the results in the form of a 
light curve and discuss possible interpretations. 

2. Observations &i Data Reduction 

Millimeter Interferometers: We used the Owens Valley Radio Observatory's millime- 
ter array and the Bcrkclcy-IUinois-Maryland Association millimeter array to observe the 
GRB 030329 in a dedicated campaign beginning 2003 March 30 through 2003 April 17. An 
additional observation was obtained on 2003 April 30. 

For the first three nights, when the redshift of the GRB was unknown, we tuned the 
local oscillator to the CS(2-1) line at a frequency of 97.98 GHz. After that, we tuned to the 
CO(J=1-0) line redshifted to z=0.1686. 

We imaged the continuum fiux using the analog correlator at OVRO which consists of 
four 1 GHz bands around the LO frequency from ±0.5-1.5 GHz and ±2.5-3.5 GHz. The array 
was in the H configuration and the projected baselines ranged from 5.5 to 84.0 kA and the 
single-sideband system temperatures ranged from 150 to 400 K. The complex instrumental 
gain was calibrated using the quasar 1159+292 every 15 minutes. The fiux of 1159+292 was 
determined using 3C273 since no planets were available for an absolute calibration; hence, 
the resulting uncertainty in the overall flux scale is ~15%. We use a flux of 1.95 Jy for 
1195+292 to cahbrate all the data. The cahbrations were done with the MMA software 
package (Scoville et al. 1993) imaged using standard routines in MIRIAD (Sault, Teuben, 
& Wright 1995). 

At BIMA we used the 800 MHz digital correlator to map the GRB. BIMA was in the 
compact (C) configuration and projected baselines ranged from 2.2-29.3 kA. The single- 
sideband temperatures ranged from 200 to 400 K. The same phase cahbrator, 1159+292 
was used to calibrate the BIMA data. The results for the continuum fiux measurements are 

summarized in Table 1. At OVRO we used the 1 GHz digital spectrometer, and the new 8 
GHz COBRA spectrometer, and at BIMA we used the 800 MHz spectrometer both to search 
for a redshifted CO spectral line in absorption; we did not detect it. 

Millimeter continuum fiux from the GRB was detected at a(J2000)=10M4™49^95, 
^(J2000) = 21°31' 17'.'30. This position is coincident with the radio and optical afterglows. 
The observing campaign is summarized in Table 1. On the first three days, we were also 
able to track the evolution of the GRB fiux within a scan because it was bright. These data 
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are listed in Table 2. The light curve is plotted in Figure 1. 

MAMBO 250 GHz observations: The millimeter continuum measurements were made 
using the 117-channel MAMBO-2 array (Kreysa et al. 1999) at the IRAM 30 m telescope 
on Pico Veleta (Spain). MAMBO-2 has a half-power spectral bandwidth of 210 and 290 
GHz with an effective frequency of 250 GHz. The beam size on the sky is 10.7 arcsec. 
The sources were observed with a single channel using the standard on-off mode with the 
telescope secondary chopping in azimuth by 32" at a rate of 2 Hz. 

Observations of GRB 030329 were done on eight different epochs between March 30 and 
April 20, with integration times on sky ranging between 5 and 17 minutes. The observing 
conditions were good, with the exception of April 2, when the sky noise was unusually strong. 
The source was observed at elevations between 56 and 74 degrees, with line of sight opacities 
at 250 GHz in the range 0.1 to 0.3. 

The pointing and focus were checked frequently on the quasar 1043+241, which is only 
3 degrees away from the GRB. Usually the pointing was found to be stable within ~ 3". For 
the absolute flux calibration a number of calibration sources were observed, mostly CW-Leo, 
which is near the target. We used a flux calibration factor of 35,000 counts per Jansky, the 
one sigma uncertainty of which we estimated to be 10%. 

The data were analyzed using the MOPSI software package. Correlated noise was sub- 
tracted from each channel using the weighted average signals from the surrounding channels. 
These data are listed in Table 3 and the light curve is plotted in Figure 1. 

3. Results 

The 100 GHz and 250 GHz millimeter light curves in Figure 1 exhibit two main features: 
a bright, constant flux density portion and a steep decline. A power-law fit to the decay gives 
= —1.98 (where F^^ cx t°) at 100 GHz and = —1.68 at 250 GHz. These values are 
in excellent agreement with that derived by Price et al. (2003) for the optical decay beyond 
~0.5 days (cto = — 1.97± 0.12) and suggests a common physical effect. 

The second feature of the light curves in Figure 1 is the bright flux density plateau. In 
the flrst week after the burst we derive a mean flux density at 100 GHz of 58 mJy. The 
emission is constant during this time expect for a two day period starting April 2 where 
there is a small (< 20%) but significant drop in the fiux density followed by a rise to its 
mean level on April 5. Likewise, the mean fiux density at 250 GHz for the first 4 epochs 
prior to the decay is 44 mJy. Small variations at 250 GHz during this time are hkely due to 
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uncertainties in the calibration. 

During the time that these miUimeter measurements were made the the optical hght 
curve had a prominent bump at At— 1.5 days with the flux increasing by a factor of two. 
There may be a second bump at 3.5 days, but this interpretation depends on whether we 
compare to a power law with a = — 2 or —1. In order to search for similar variations at 100 
GHz we subdivided each of our observing sessions into 2-hour intervals. The results of this 
analysis are summarized in Table 2. With the exception of a possible increase by about 10 
mJy (15 — 20%) at t ^ 0.7 d and t ^ 2.8 d, neither of which corresponds to a change in the 
optical brightness, we find no significant variations. 

4. Discussion 

Two explanations have been offered as to why GRB 030329 appears to be under- energetic 
compared to other GRBs. Granot ct al. (2003) have explained the observed steep decay in 
the optical light curves of GRB 030329 at t = 0.55 days (Price et al. 2003) in terms of a 
standard model of a jet expanding into a constant density medium (Sari, Piran & Halpcrn 
1999). Note that such steep temporal breaks Act > 1 are not expected in spherical outflows 
in general, nor are they expected in coUimated outflows expanding into a wind-blown medium 
(Sari et al. 1998; Kumar & Panaitescu 2000). In their model the fluctuations in the optical 
light curve (§3) originate from "refreshed shocks", in which slower moving ejecta catch up 
with the main shock and re-energize it. In this picture energy injection from the central 
engine is not instantaneous (as commonly assumed) but episodic, with most of the energy 
being carried by ejecta with the lower Lorentz factors. 

An unavoidable consequence of having the GRB central engine inject a range of Lorentz 
factors r over time, instead of a single high Lorentz factor shell, is that each newly arrived 
shell at the shock front modifies the afterglow spectrum, producing a short-lived reverse 
shock, and shifts the combined spectrum to lower frequencies (Rees & Meszaros 1998; Kumar 
& Piran 2000). According to the more detailed calculations of Sari & Meszaros (2000), which 
take into account synchrotron self-absorption, there should be a significant enhancement of 
the flux density at millimeter and submillimeter wavelengths. Short-lived reverse shocks 
increase the peak flux density by F and shift the peak to lower frequencies by a factor F"^. 
For the Lorentz factors given by Granot et al. (2003) we would have expected to see order of 
magnitude variations in the millimeter flux decaying on timescales of a day or less. This key 
prediction is contrary to what was seen (§3) in our millimeter data and thus our observations 
rule out the refreshed shock model. 
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Berger et al. (2003) proposed a two component jet model based on the existence of two 
different jets breaks, one in the optical light curves at 0.55 days and another best seen in the 
radio light curves at 9.8 days. Single-jet fits to the radio data could not explain the evolution 
of the optical emission before 1.5 days, especially the sharp break at 0.55 days (Price ct 
al. 2003). The millimeter observations presented here are crucial in this respect since they 
define the peak of the synchrotron spectrum and hence the overall normalization. Lacking 
such data, earlier claims of a two component outflow for GRB 991216 (Prail et al. 2000) are 
less secure. In this case Berger et al. (2003) define a narrow angle jet with t^Aj — 0.55 
days and 9naj — 0.09 rad which is responsible for the early afterglow, and a wide angle jet 
twAj — 9.8 days and 9waj — 0.3 rad which carries the bulk of the energy in the outflow and 
dominates the optical and radio emission after ~ 1.5 days. It is likely that the plateau, seen 
in the millimeter light curves (§3) during the first week, is a combination of the falling fiux 
density (as t~^/^) from the NAJ, and a rising fiux density (as t^/^) from the WAJ (Sari et 
al. 1999). 

Since in the two-component model the fiux increase at 1.5 d is due to the rise of the 
wide-jet component, which then evolves in the usual fashion as F,y oc t~^, the subsequent 
possible fiuctuations discussed by Granot et al. (2003) and in §3 are modest and no longer 
require refreshed shocks but instead could be explained by variations in the circumburst 
density. 

In summary, the millimeter observations presented here have been used to distinguish 
between two equally compelling models for GRB 030329 and its afterglow. The absence of 
bright, short-lived millimeter emission, coincident with "bumps" in the optical light curve 
between 1 and 7 days after the burst, was used to show that the GRB central engine was 
not actively injecting energy on this timescale. Instead, the millimeter data support the 
proposed two-component jet. It is possible that the true structure of GRB outfiows are 
considerably more complicated than the simple picture presented here (Zhang & Meszaros 
2002; Rossi, Lazzati & Rees 2002; Perna, Sari & Frail 2003). Models of relativistic jets 
propagating out through the stellar progenitor show that there may exist a large range 
in Lorentz factors in the outfiow (MacFadyen, Woosley & Heger 2001; Zhang, Woosley 
& MacFadyen 2003), which decrease away from the rotation axis as the degree of baryon 
entrainment increases. Nonetheless, within the limits of the current data for GRB 030329, 
there is evidence for jet structure with at least two distinct components, with the wider of 
the two carrying the bulk of the energy. This last point is worth emphasizing since events 
hke GRB 030329/SN 2003dh and GRB 980425/SN 1998bw (Kulkarni et al. 1998) have shown 
that the group of sub-energetic bursts may simply be an artifact of limited observations. True 
calorimetry of GRBs must account for material at low Lorentz factors (Berger et al. 2003), 
which typically are brightest at radio and optical wavelengths. 
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Fig. 1. — Millimeter flux as a function of time on a logarithmic scale. The stars indicate 100 
GHz (3mm) data. Filled stars represent OVRO observations, and the open stars represent 
BIMA observations. Filled circles represent the 250 GHz, IRAM 30-m MAMBO observa- 
tions. On two days (2003 April 10 and April 11), both the BIMA and OVRO observatories 
observed the GRB. The measurements are consistent with each other. The inset shows 
the measurements from the first two OVRO epochs which have been subdivided into 2-hr 
intervals to look for flux density variations. 
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Tabic 1. Observations 
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Table 2 


Early Time Evolution 
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Table 3. 250 GHz MAMBO Observations 
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